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We establish the model-independent spin structure of the matrix elements for the near-threshold 
scalar meson production in pp— and np-collisions, when the final particles are emitted in S-state. 
Polarization phenomena are derived in a general form. The properties of the t— channel dynamics, 
which is based on different meson exchanges, are studied in terms of the s— channel parametrization 
of the matrix element. The prediction of a 'realistic' model, based on n + cr-exchanges are also 
£S| ' presented. 
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The strong and electromagnetic decays of the scalar mesons S (S = a, /o, and do) are recently object of large 
theoretical and experimental interest. The structure of these mesons is not yet fully understood and the coupling 
constants of their decays not determined (see |l],| ,||j|f]Uf7[H and refs. herein). As an example, the electromagnetic 
constants g p<T1 and g W(T7 are very important for the solution of different problems in in hadron electrodynamics. In 
the near threshold region, the constant g pai drives the cr-contribution to the differential cross section of the process 
7+p —>■ p + p + p a ||[l0 f 1 1 • The coupling constant g^a~i is important for the estimation of the effects of the 



p meson exchange currents for the deuteron electromagnetic form factors [ [l 2|Jl3|Jl4f| , particularly in the region of large 

momentum transfer. Both these constants play a major role in the interpretation (f5) of the HERMES effect [ JL6| , 
concerning the anomalous behavior of the electroproduction cross section on nuclei at low Q 2 - where the cross 
O^N section is enhanced for longitudinally polarized virtual photons and depleted for transversally polarized photons. 
Moreover, the constants g pcr -y and g^a^ enter in the interpretation of different radiative decays of vector mesons, like 
p°(uj) -> 7r°7r°7 and p°{u) -> 77+71-7 0. 

Two-photon decays of scalar mesons, S — > 27, which are important for the estimation of the corresponding t— channel 
' contributions to the amplitudes of real and virtual Compton scattering on nuclcons |T^,|I^,^,^f) , could be estimated 
1 , on the basis of the gva-y coupling constants, in the VDM approach. 

In the framework of the effective Lagrangian approach, the calculation of the cross section for scalar meson photo- 
and electro-production on nucleons requires the knowledge of the gva-y and gjvws-coupling constants. The strong 
coupling constants gNNS |2^ , |23[| enter in various calculations in hadron dynamics, not only with respect to the NN- 
potential, but also for different observables for processes like tt + N -> S + N N + N — > S + d p5|p^,|27|,p8|| , 

and N + N -► S + N + N ||. 

The feasibility of the experimental study of scalar meson production depends essentially on the nature and on the 
rate of the decay. If, for example, the decay fo — > KK or ao — > KK dominate, these mesons could be observed in the 
KK effective mass distribution, close to the kaon mass, as a resonant contribution. Such effects have been observed 
in 7riV-collisions P, p9tp^0l|3l|j3^ , |33| , ^| . The interpretation of the data needs an adequate theoretical approach for the 
process 7r + N — > N It was shown that the reaction p + p — > p + p + K + K~ in the kinematical conditions of the 

DISTO |35| or at COSY the presence of the f signal is hindered by a large background. Therefore, in principle, the 
process n + p — * n + p + K + K~ could be more favorable, because the cross section is one order higher, whereas the 
background is comparable with respect to pp-collisions. 

Note that the experimental study of if + A"~-production in the process p + p^p + p + K + + K~ , at an energy 
excess Q=17 MeV over threshold j36fl, gives the following value: o{pp — » vvfo) — (1.84 ± 0.29^035) nb, including the 
statistical and systematic errors. The extension of this study is foreseen [|37| , p8|j39| l . 

In this paper we derive the most general and model independent properties for the processes of scalar meson 
production in NN— collisions, in the threshold region, where the theoretical analysis is essentially simplified. The 
spin structure of the corresponding matrix elements contains a set of sixteen independent amplitudes in the general 
case, eight amplitudes for coplanar kinematics and only three independent amplitudes in the threshold region. 
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This paper is organized as follows. In Section II, using the selection rules with respect to Pauli principle, P-parity 
and total angular momentum, we establish the spin structure of the matrix elements for the processes p+p — > p+p+S° 
and n + p — > n + p + S° and analyze the polarization phenomena for these processes in a model independent way. In 
Section III we transform (using the two-component Fierz transformation) different t— channel contributions, which are 
described by a definite set of Feynman diagrams, to the universal s-channcl parametrization and find the expression 
for the corresponding partial threshold amplitudes. In Section IV we discuss the predictions of a realistic model based 
on a + 7r-exchanges. A short discussion of final state interactions (FSI) is done in Section V. In the Conclusions we 
summarize the results obtained. 



II. SPIN STRUCTURE OF THE THRESHOLD MATRIX ELEMENT AND POLARIZATION 

OBSERVABLES 

The spin structure of the matrix elements for scalar meson production in NN— collisions, N + N^N + N + S°, 
in the threshold region is determined by the selection rules with respect to P-parity, total angular momentum and by 
the Pauli principle. The threshold region, where all final particles mN + N^N + N + S° are produced in relative 
S-state, can be rigorously described by a formalism based on the two-component nucleon spinor parametrization of the 
corresponding matrix element. This formalism can be built in model independent way MOf, and has been previously 
successfully applied to vector Q] , pseudoscalar [Q , and strange particle production in nucleon- nucleon collisions. 

Let us explicitly derive the matrix element for the processes p + p — > p + p + S° , where S° denotes a neutral scalar 
meson, S° = a, fo or ajj. Due to the Pauli principle, at the reaction threshold, only quantum numbers j p = + are 
allowed, where j is the total angular momentum and P is the P-parity of the colliding protons. Only a single partial 
transition can take place, corresponding to: 

Si = 0, I = -» f = 0+ -> S f = 0, (1) 

where Si (Sf) is the total spin of the initial (final) protons and I is the angular orbital momentum of the colliding 
protons. 

The spin structure of the threshold matrix element for the transition in the CMS of the considered reaction, 
can be parametrized in the following general form: 

M pp = g(xlvyxl)(X2VyXi), (2) 

where xi an d X2 (X3 an d Xi) are the two-component spinors of the initial (final) protons; g is the threshold partial 
amplitude, describing the singlet-singlet transition in the pp-system (with scalar meson production in the S'-state), 
which is generally a complex function of three independent energies: W (the total invariant energy of the colliding 
particles), E\ and E 2 (the energies of the scattered protons). So, all the dynamics of the considered process is 
included in the amplitude g, but the exact form (^|) of the spin structure of the matrix element results from a 
generalized quantum mechanical kinematics. 

The Pauli matrix a y , in the parametrization (^) , insures the correct transformation properties of the corresponding 
two-component spinor products, relative to rotation. 

The presence of a single amplitude in (|^) implies that the spin directions of all the protons are fixed, with definite 
relative angles and fixed modules. In other words, all polarization phenomena for p + p—>-p + p + S at threshold can 
be predicted without knowing the amplitude <?, i.e. in model independent way. Moreover, these polarization effects 
are the same for any scalar meson: a, fo or a{j. However, the absolute value of the cross section, which depends on 
\g 2 \, is different for different scalar mesons. 

All T-odd polarization observables (i.e. one-spin and three-spin polarization correlations) are identically zero. The 
dependence of the cross section on the polarizations Pi and P2 of the colliding protons can be written as: 

dio \doj J \ / 

(typical for the singlet pp— interaction), where doj is the phase space volume element for the three-particles final state 
and (da/duS)o is the differential cross section with all unpolarized protons in initial and final states. 

All polarization transfer coefficients, characterizing the dependence of the polarization of any final proton from the 
polarization of the initial proton, vanish also. 

The situation is very different in case of np— collisions, n + p—>-n+p + S , in the threshold region, where there 
arc three allowed transitions: 
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The corresponding matrix element can be written 
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the following 


form: 
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-fj2 



Xl(<*a ~ k a O ■ k)<Jyxl {X2<?y{Oa ~ • fc)*! 



(3) 



where k is the unit vector along the three-momentum of the initial neutron beam, a = (a x , a y ,a z ) is the standard set 
of Pauli matrices, and g\ — gs are the partial amplitudes for the rip-interaction. Due to the isotopic invariance of the 
strong interaction the following relation holds: 



9i 



7 9 , 



(4) 



i.e. the triplet-triplet amplitudes g% and g% are present only in np-collisions. Therefore polarization phenomena in 
n+p—>n+p + S° collisions are more complicated then in pp-collisions. There is no universality here, polarization 
phenomena are different for the different scalar mesons. There is, however one general feature: all one-spin observables 
vanish, for n + p — > n + p + S°, in particular the T-odd polarization of the final nucleons, emitted in the collisions of 
unpolarized nucleons and the analyzing powers in n + p- or n + ^-collisions. 

The dependence of the differential cross section for n + p-collisions can be written in the following general form, 
which is correct near threshold: 



da 



(5) 



where the real spin correlation coefficients A\ and Ai are determined by the following formulas (in terms of the partial 
threshold amplitudes <?i, i = 1 — 3): 



A 2 



da 
du) 
da 
du) 



= - l.9i I + |S3| , 

= 2(| 52 | 2 -| 53 | 2 ). (6) 
The amplitudes gi are normalized in such way that: 

M 2 + 2| ff2 | 2 + | 53 | 2 , (7) 
where (da / 'du>)o is the differential cross section with unpolarized particles. So, from Eqs. (^|) and (Q) one can find: 



da 
du> 



4| 91 | 2 = (1 - 3Ai - A 2 ) 



da 

du 



4|//2|- = n ■■!■■ A ^+ A ^\% 



4|//,f = (] +.4^.4,11^ , 



dui J 

which shows that the moduli of the three threshold amplitudes for n + p — > n H 
measurement of Ai, A2, and (da/duj) . Taking into account the isotopic relation 



(8) 



(- S° can be determined by the 
one obtains: 



K 



(da np /duj) 
{da pp /duj) 



= 1-A 1 -A 2 
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The dependence of the polarization P 3 of the final neutron on the polarization Pi of the initial neutron can be 
parametrized as: 

P 3 = P1P1 + p 2 k(k ■ A), (9) 
where p\ and p 2 are the spin transfer coefficients, which can be expressed in terms of the partial amplitudes gf. 

2neg 2 ( 9l -g 3 y 2\g 2 \ 2 - 2Ue [ gl g* - g 2 ( 9l - g 3 )*} 

\ 9l \2 + W + \g 3 \^ P2 l3i| 2 + 2| 52 | 2 + |.g 3 ' 2 



The dependence of the polarization P4 of the final proton on the polarization Pi of the initial neutron can be also be 
parametrized in terms of two real coefficients: 

P4 = P3P1 + P4k{k ■ Px) (11) 

with the following formulas for p 3 and p^: 

2Keg 2 (gi+g 3 y _ 2\g 2 \ 2 + 2Ke [ gi g 3 * - g2 ( gl + .93)*] , , 

P3 | g i| 2 + 2| 52 p + | 53 | 2 ' Pi l5i| 2 + 2| 52 | 2 + M 2 " 1 ) 

Comparing Eqs. (JsJ) , (]?]), (|To|), and (|lj), one can find the following relation between these polarization observables: 

Pi+P2+p 3 +P4 = l + A 1 +A 2 , (13) 

which has to be verified by any model describing this process. 

Note that the isovector ao-meson can be produced in the following processes of the ./ViV-interaction: 

p + p^ p + p + Q,Q, 

p + p—> n+ p + Oq , 
n + p—> p + p + a,Q , 
n + p—> n + p + a,Q. 

It is possible to show that the production of the charged aj-mesons -in the near threshold regime is characterized by 
a single amplitude, with the spin structure as in Eq. (|^), where the following relations hold: 

g(pp -> npa%) = -g{np -> ppa^ ) = -^9i n P PP a o)- 

Therefore the previous statements on polarization phenomena for p + p—>p + p+S° apply also to the production of 
the charged -mesons. 



III. THE DYNAMICS FOR THE t-CHANNEL 



The standard dynamics for different processes of meson production in NN-collisions, N + N — > N + N + V , 
N + N — > N + N + P, and N + N — * N + N + S° , (where P, S, and V denote pseudoscalar, scalar and vector mesons, 
respectively) is based on the consideration of mesonic exchanges in t— channel, such as tt, 77, cr, ao, p, to, etc. In such 
approach one has to know the meson-nucleon coupling constants, PNN, SNN, and VNN and the amplitudes of 
different subprocesses, such as : 

P* + N + S, 

S* + N + S, (14) 
V* + N + S, 

where the index '*' denotes virtual mesons - with space-like four-momenta. In principle some information exist on 
the coupling constants, but the threshold amplitudes for the processes ([hf ) are poorly known, in particular for those 
processes that can not be experimentally studied. Therefore model calculations have to be done, in order to find these 
amplitudes. In case of complex amplitudes not only the absolute values are important, but also the relative signs and 
phases. 

In this section we analyze different t— channel exchanges for the processes of scalar meson production in AW- 
collisions, and give, when it is possible, model independent predictions. In other words, we will find expressions for 
the polarization phenomena mn + p —> n+p + S°, which depend only on the quantum numbers, spin and parity J p , 
and isospin T, of the exchanged meson, but not on the corresponding coupling constants, hadronic form factors and 
threshold amplitudes. Let us consider the spin structure of the matrix element (0) for different exchanged particles. 
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A. Scalar exchange, J p — + , 1=0 

Taking into account the identity of initial and final protons in p+p — * p+p+S°, the a-exchange here is characterized 
by four different t— channel diagrams (Fig. Q). Each diagram has a different spin structure, with different order of the 
two-component spinors Xii i = 1 — 4. Only the sum of all these diagrams, which satisfies the Pauli principle for the 
initial and final protons, generates the correct spin structure (||), with a single amplitude g. In threshold conditions 
the spin propagator for the considered diagrams is written as: 

I I 



t — m 2 a mm s + w? a 

where m, m a , and m s are the masses of the nucleon, of the virtual cr-meson, and of the produced scalar meson, 
correspondingly. Such equality is correct at the reaction threshold, and it follows from the assumption of 5*-wave 
production of final particles. The difference in these propagators, which appears far from threshold, generates P- and 
higher waves of produced particles. 

The er-exchange for n + p — ► n + p + S° is described by two diagrams, only (Fig. ^|). Taking into account the 
isotopic invariance of the strong interaction, one can write the following threshold matrix element for the process 
n + p — > n + p + S°, corresponding to cr-exchange: 

M =Mla+M 2 a = 2Mi <7 , 

M la = ^LMA{a*p -> Sp)(xll X 2)(xUxi), (15) 

where g a NN is the aNN coupling constant, X is the unit 2x2 matrix, A(a*p — > Sp) is the partial amplitude for the 
subprocess a* + p ^ S + p - in the S— state, J\f = 2m(Ei + m) = m(Am + m s ) is the normalization factor, related to 
the transformation from four-component Dirac spinors to two-component Pauli spinors. 

In order to find the partial amplitudes gi, corresponding to cr-exchange, it is necessary to apply the Fierz- 
transformation, in its two-component form, to the spinor construction in Eq. (|f5|): 

(xUx2){xUxi) = -\[x\^yX%){x^yXl) + ^(xW<Tyxi){X2<Jy<TaXi)- 

Comparing this with the general parametrization of the spin structure for the threshold matrix element of the process 
n+p—>n + p + S a , one can find the following expressions for the particle amplitudes c^: 

gi a (np -> npS°) = -A a , 

g2a{np -> npS°) = A a , (16) 
93a{np -> npS Q ) = A a , 
g a (pp - ppS°) = -2A a , 



where 



A a = NA(ap — > pS°)- gcrNN 



t 



Eqs. (|l6|) allow to predict all polarization phenomena in n + p — * n + p + S°, independently on the details of the 
considered model, such as the value of the constant g a NN and the amplitude A(ap — ■> pS ). 
One can find: 

Pier = 1, A la = A 2a = p 2a = P3a = Pia = 0, (17) 



Ua = ^("P -> n P S °) = 2 (18) 
a <y{pp —> ppS°) 

In (|lj) we took into account the identity of the produced protons in the reaction p + p^p + p + S°. 
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B. ^-exchange, J F = , 1=0 

As in the case of er-exchange, the ^-exchange generates four different diagrams for pp-interaction and two diagrams 
for np-interaction. Let us consider, therefore, the simplest case of n+p—>-n+p + S^-processes, the matrix element 
of which can be written as: 

M = Mir, + Mir, = 2Mir„ 

Mir, = -A v (xla-kx2)(xlv-kxi,) (19) 



hi + m t — Tnfi 



JVk 

E + n 

where g n NN is the coupling constant for the rjNN-vertex, A(t]*N — > NS°) is the threshold amplitude for 77 + N — > 

iV+S 10 , which describes the partial transition: £„ = 1 —* j p = 1/2+. The kinematical factor — = ./ (at 

' £ + m V 4m + m s v 

threshold) results from the transformation of the pseudoscalar vertex to its two-component equivalent x^cr • k\- 
After applying the Fierz transformation to the spin structure ( |l9| ) to the standard form, Eq. (^), one can find the 
following formulas for the partial amplitudes gi V , i = 1 — 3, corresponding to 77-exchange: 

— > npS°) = -A v , 

g 2v (np^n P S°) = -A r „ (20) 
93r,{np -> npS°) = A,, 
9r,{pp -> PPS") = -2A V , 

and similarly, to cr-exchange, one can predict all polarization observables: 

Air, = Air, = P37] = Pir, = 0, Plr, = -1, P2r, = 2,11 = 2. (21) 

Comparing Eq. ( |l7| ) and Eq.(^), one can see that the polarization transfer coefficients, characterizing a change of 
neutron polarization, must be very sensitive to the quantum number of the isoscalar exchange in t— channel, other 
polarization observables, such as A\, A2, P3 and p^ vanish for both exchanges J p = + and J F = 0~. 

C. Pion exchange, J p = 0~, 1=1 

Any isovector exchange for the process n+p — > n + p + S° is characterized by a set of four Feynman diagrams Fig. 
(|J) with exchange of neutral and charged pion. 

Considering the contribution of all these diagrams, and taking into account the isotopic invariance predictions for 
the two vertices of the considered diagrams: 



1 1 



one can find (after applying the Fierz transformation) the following expressions for the partial amplitudes g^, i = 1 — 3, 
corresponding to 7r-exchange: 

giir{np — > npS°) = A v , 

g 27C (np -> npS°) = -3^, (22) 
5377 (np -> npS°) = 3A V , 
gAPP -> PP S °) = 2 An, 



where A n = AfJ — — 9* NN A{n*N -> NS). In principle, the amplitude A(tt*N -> NS) for threshold S- 

y 4m + m s t — m£ 

production in 7riV-interaction can be experimentally measured through the corresponding differential cross section, 
but this method will not determine its sign. 



G 



However numerical values of the polarization phenomena for the n + p — > n + p + S° processes can be predicted 
exactly, without information on A(ir*N — > NS): 

A l7T = 2/7, A 27r = 0, p l7r = -3/7, p 27[ = 6/7, p 37r = -6/7, p A „ = 12/7, (23) 

and they are very different from the case of isoscalar exchange. 

Another interesting result concerns the large difference in cross section, for pp- and rip-interaction: TZ n = 14, i.e. 
for scalar meson production (in NN— interaction) we have large isotopic dependence. 



D. Scalar exchange, J F = + , 1=1 

In this case, both processes p + p^p + p + S° and n + p^n+p + S° is characterized by a set of four Feynman 
diagrams, the neutral ao-exchange for pp— collisions and the neutral+charged ao-exchanges for np-collisions. 
The partial amplitudes gi ai i = 1 — 3, corresponding to ao-exchange are determined by the following formulas: 

gia(np -> npS a ) = A a , 

92a{np -> npS°) = 3A a , (24) 
33a (np -> npS°) = 3A a , 
9a (pp -> ppS°) = 2A a , 

and A a = Af 9aNN A(a*N -> NS). The threshold amplitude A(a*N -> AT 5) for the exotic subprocess a* + N -> 
t - 

A + .S has to be determined in framework of a model. 

From these formulas, independently from the values of the coupling constant g a NN and the threshold amplitude 
A(a*N — > iVS), one can predict the following values for the polarization observables: 

Ai a = 2/7, A 2a = 0, pi a = 3/7, p 2a = 0, p 3 a = 6/7, p 4a = 0, ^ a = 14. (25) 

They are different from the case of isoscalar exchange and from (]25|). Due to isovector exchange, the isotopic ratio 1Z 
is large also in this case. 



E. Correlation of spin and isospin structures for the threshold regime 

Combining the previously derived contributions, one can find for the threshold amplitudes <?.; the following expres- 
sions: 

gi (np -> n P S°) = -A a + A a - A v + A v , 

g 2 {np^npS°)=A a + '3A a -A v -3Ai T , (26) 
g 3 (np -> npS°) = A a + 3A a + A v + 3A V , 
g(pp — > ppS a ) — 2gi(np — » npS°). 

Let us note that the S-channel parametrization of the spin structure for the threshold matrix elements allows 
(after applying the Fierz transformation) to unify different t— channel contributions- with different J p - in universal 
and transparent form, which is well adapted to the analysis of the sensitivity of the polarization phenomena, in 
n + p — > n + p + S°, to the quantum numbers of the t— channel meson. Such unification allows to simplify all 
calculations: we can express the matrix element for any exchange in terms of three amplitudes only, whereas, for 
example, a model with r\ + 7r + a + a contains 2 + 4 + 2 + 4= 12 different Feynman diagrams - with different spin 
structures. All these twelve contributions to the total matrix elements, in general, interfere, so it is in principle 
necessary to calculate 12 x 12 = 144 terms, instead than 3x3 = 9 terms, for the gi parametrization of the matrix 
elements. Moreover, these nine terms are the same for any model, whereas, for example, adding vector exchanges 
will increase essentially the number of t— contributions. Another advantage of the t — > s Fierz transformation is the 
explicit dependence of the observables on a definite combination of coupling constants, hadronic form factors and 
elementary amplitudes. This helps in finding out which contributions play the most important role in the t— channel 
dynamics, and gives a feedback on the coupling constants by comparison with the experimental data. 

Another important property of the partial amplitudes gi for np— processes, Eq. (|2^), concerns a strong correlation 
of the spin and isospin structure of the matrix element in the threshold region. The amplitudes g 2 (np — > npS°) and 
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gs{np — > npS°), describing the np-interaction in the isotopic singlet state, are in general different, even in case of 
a definite isospin in t-channel, if the P-parity of the ^-channel is not fixed. However, for definite J^-exchanges, we 
have: 

g 2 (np -> n P S°) = g 3 (np -► npS°), if J p = 0+, X = and 1=1, 

# 2 (np -> np5°) = -# 3 (np -> if J p = 0~, X = and 1=1. (27) 

The relations ( p7|) hold also for any combination of isoscalar and isovector exchanges. 

The total cross section for n + p — > n + p + S° is the incoherent sum of isosinglet and isotriplet contributions: 

a{np — > npS°) — ao{np — > npS°) + a\{np — > npS°), 

a Q (np -> npS°) = 2| ff2 (np -» npS°)| 2 + | ff3 (np -» ripS )] 2 , (28) 
(Ji (np -> npS* ) = -> r7,pS°)| 2 = -cr(pp ppS* ). 



IV. A POSSIBLE MODEL FOR N + N^N + N + S a 

We analyze here in detail a model for N + N — > 7V + N + S°, based on cr + 7r-exchange. In order to justify such model, 
let us mention that the 77-contribution can be neglected, in (E6|), due to the fact that the coupling constant g V NN is 
presently poorly known, being in the range l-j-7 |4^^5| , ^ , [47| | . Concerning the ao-exchange, both the ingredients of 
such contribution, the coupling constant g a NN and the amplitude A(aN — > NS°) are presently not known. The mass 
of the ao-meson is larger in comparison with the pion mass. In any case, the cx + 7r-model can be considered a realistic 
starting point for the analysis, containing the exchange of mesons with different P-parities and isospin. As a result, 
all three partial amplitudes are different: 

g\(np — > npS°) = ^(1 - r), 

g 2 {np^npS°)=A v (-3 + r), (29) 
53 (np -> npS°) = Ax(3 + r), 
g(pp^ P pS°)= 2Ar(l-r), 

where r is the ratio of the corresponding contributions: 

gaNN (t-ml\ A(a*N^NS° 



r 



g-xNN \t-mlJ A(ir*N -► NS°) 




One can see that all the physics of this model (with eight different Feynman diagrams) is contained in a single 
parameter r, which is basically the ratio of the coupling constants and the elementary amplitudes. In the general 
case, the ratio r is a complex parameter, which depends on the excitation energy of the produced NNS-system. 
Therefore the polarization phenomena for the process n + p — » n + p + S° and the ratio of cross section for np- and 
pp-collisions, can be expressed in terms of two parameters |r| 2 and IZer: 

a(np^npS°) 7-2^er+|r| 2 / |r| 2 - 2Ker \ 

a(pp^ppS°) l-2Ker+\r\ 2 \ 1 - 2Ker + \r\ 2 J ' 1 ' 

The coefficients Ai and A2, which characterize the polarized np-collisions, can be written as follows: 

A — 9 l+7£er ■ _ 3TZer 

1 ~ 7-2fter+|r| 2 ' 2 ~ 7-2fter + |r| 2 ' ( j 

The ratio r can be found in framework of a model for the elementary subprocesses it* (a*) + N — > N + S° , under 
several assumptions concerning coupling constants, cut-off parameters, form factors... To avoid the uncertainties 
related to these choices, we assume that r is real. This is the case, for example, of the effective Lagrangian approach, 
which gives real amplitudes and, therefore, real values for the ratio r. Another possibility is the saturation of the 
A(a*N — ► NS°) and A(n*N — > NS°) amplitudes by a single nucleon resonance with J p = l/2 + . Such mechanism 
produce complex amplitudes with zero relative phase, and r is a real parameter, too. 
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In this case we can predict the r-dependence of 1Z, Ai, and A2 (Figs. 4-6). The ratio R of the corresponding cross 
sections can be very large in the region r ~ 1 (Fig. 4). The experimental determination of R will allow to find two 
solutions for r: 

The function A\ (r) has two extrema: the maximum at r = — 1 + \/lQ, where A" lax — (2 + •\/T0)/v / 6 and a minimum 
at r = -1 - \/l0, where Af in = (2 - v 7 !?))/^, see Fig. 5. 

The function A^ir) has a very similar behavior, with a the maximum at r = +v7, where A™ ax = (1 + \/7)/4 and 
a minimum at 7- = — \/7, where A™ 1 ™ = 12— v7)/4, see Fig. 6. 

V. COMMENTS ON THE NN FINAL STATE INTERACTION (FSI) 

In order to illustrate possible effects of FSI for the considered reaction, let us make an oversimplified estimation, 
using only the TV TV-scattering length approximation, which will result in an upper limit for the NN-FSI. Let us 
consider, for illustration, the case of pion exchange. One finds the following correction for some observables which 
have been previously discussed: 

• The ratio of np and pp-cross sections: 



1 


~a s {np)~ 


2 27 


'a t (np)~ 


2 


a s{PP)_ 


+ y 


a s(PP)_ 



• Polarization effects for np-collisions: 

nn 27 + R 2 st ' nn 27 + R 2 st ' 

where R s t = [a s (np)/at{np)], a s (NN) is the singlet iViV-scattering length and at(np) is the triplet np-scattering 
length (for S-state). 

Using the following values for the scattering lengths Jig ): et s (np)=-23.768 fm, a t (np)=5A2A fm, and a s (pp) =-7.8098 
fm, one can find: ijW ~ 11 (instead of 14), C$ ~= -0.22 (instead of 0.29), and D ( nJ --0.96 (instead of -0.43), i.e. 
the final state np-interaction can change the sign of Cir? , but does not affect much the ratio R^ . 

VI. CONCLUSIONS 

Let us summarize the main results concerning the production of scalar mesons in nucleon-nucleon collisions near 
threshold. 

• We parametrize the spin structure of the threshold matrix element for the processes p+p — > p+p + S°, where S° 
is a neutral scalar meson, S = <r, /o, or ao, in terms of a single spin structure, corresponding to a singlet-singlet 
transition in the pp— system. The same structure describes the production of the charged a^-mesons in the 
processes: p + p^p + n + a,Q and n+p^p + p + a^. 

• The process n+p^p + p + S° is characterized by a more complicated spin structure, where the matrix element 
contains, in the general case, three contributions with different spin structures. One contribution coincides with 
the Ai(pp — > ppS°) matrix element -due to isotopic invariance, and two additional contributions describe the 
triplet-triplet transition in the np-system, which are forbidden by the Pauli principle, in case of the pp— system. 

• The suggested model-independent parametrization of the spin structure, which is based on the general symmetry 
properties of the strong interaction ( such as P-invariance, isotopic invariance, conservation of angular momentum 
and Pauli principle) allows to analyze polarization effects for pp- and np— collisions in a transparent way. These 
effects are very peculiar for threshold conditions, where all final particles are emitted in relative ^-states. 
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• The standard t— channel dynamics for the process N + N — > N + N + S° is generated by different meson 
exchanges with J F = + and J p = CP and isotopic spin 1=0 and 1=1. The dependence of the polarization 
phenomena in n + p — > n + p + S° on the quantum number of the exchanged mesons can be described in terms of 
the corresponding S-channel partial amplitudes, with the help of the Fierz transformation in its two-component 
form. 

• As an example of a 'realistic' model in the threshold region, we considered the 7r + er-exchange. An attractive 
property of such model is that the polarization phenomena in np-collisions depend on a single complex parameter 
r, which characterizes the relative role of a- and 7r-exchange. The polarization observables, which do not vanish, 
show a large sensitivity to r. 

• The processes p + p^p + p + S° and n + p — > n + p + S° show large isotopic effects, i.e. S-state production of 
scalar mesons is larger in rip-collisions than in pp-collisions, and polarization phenomena are different for these 
processes. 

Note that these results do not depend on the structure of the fo and ao-mesons (gg-states, four-quark state, K K- 
molecule etc.) and on the mechanisms of their decays. 
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FIG. 1. Feynman diagrams for cr-exchange, for p + p^p+p + S°, where p(xi), i= 1-4, means that the corresponding proton 
(in initial or final state) is described by the 2-component spinor \i- 
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FIG. 2. Feynman diagrams for cr-exchange for n + p —> n + p + S°. Notations as in Fig. 1. 
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FIG. 3. Feynman diagrams for 7r-exchange for n + p —> n + p + S°. Notations as in Fig. 1. 
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FIG. 4. Dependence of R on r, see Eq. 30 
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FIG. 5. Dependence of Ai on r, see Eq. (31) 
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FIG. 6. Dependence of A2 on r, see Eq. (31). 
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